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Abstract. We remark that the power diagrams from computer science are the spines of 
amoebas in algebraic geometry, or the hypersurfaces in tropical geometry. Our concept 
of a Morse poset generalizes to power diagrams. We show that there exists a discrete 
Morse function on the coherent triangulation, dual to the power diagram, such that its 
critical set equals the Morse poset of the power diagram. In the final section we use 
Maslov dequantization to compute the medial axis. 



1. Introduction 

Power diagrams were introduced by Aurenhammer in |Aur87j , as a generalization of the 
more commonly known Voronoi diagrams. They recently surfaced in algebraic geometry, 
as the spines of amoebas, see |PR04j . Another appearance of the power diagram is what 
other authors call a tropical hypersurface, see [RGSTOSj. They also appear when string 
theorists use toric geometry, see for instance [DFG02j . 

The treatment they receive in [PRQ4j is different from the one in |Aur87j . The essential 
tool that Passare and Rullgard use is the Legendre transform, a tool that first seems to 
have been introduced in convex analysis in |Man39j . The essential tool that Aurenhammer 
uses is that the power diagram of some points in M" is a projection of a convex hull in 
]^n+i_ spell out the relation between the classical theory of Voronoi diagrams and 
that of the tropical hypersurface. Note that what we call here a Delaunay triangulation 
is called a coherent triangulation in chapter 7 of |GKZ94j . Such triangulations are called 
regular in |Zie95j . 

After that we will explore the relationship between Voronoi diagrams and power dia- 
grams. In particular given a power diagram T of A points {Pi, ■ ■ ■ , P/v} in we will 
construct a Voronoi diagram T of A points in M"+^ a hyperplane section of which is the 
diagram T. Thus power diagrams are affine Voronoi diagrams in a very precise sense. 

Then we focus on higher order structures. With a Voronoi diagram there come the 
k-th order Voronoi diagrams and the Delaunay triangulation. The k-th order power 
diagram is again a power diagram, see |AI87j . and to it there is associated the k-th order 
coherent triangulation. The union of all the power diagrams is a hyperplane arrangement. 
We conclude that not every picture that looks like a power diagram is actually a power 
diagram. 
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The main new point of view we present in this article is the relation between tropical 
hypersurfaces and discrete Morse theory introduced by Forman in |For98j . This seems to 
have gone unnoticed before. The triangulations are polyhedral complexes, and the affine 
linear functions that define the power diagram give rise to discrete Morse function on 
the dual simplicial complex. We apply this third way of looking at power diagrams to a 
notion we introduced before, the Morse poset. Morse posets are seen to be determining a 
discrete vectorfield. 

We apply our newly found knowledge to the problem of enumerating different Morse 
posets. 

In the final section we show another application of tropical geometry to computational 
geometry. We use Maslov dequantization to determine the medial axis. 

Many of our results can in some explicit form or another be found in the literature. We 
have gone to great lengths to include appropriate references, but in all likelihood we are 
not complete. It has been remarked by other authors also that many results in tropical 
geometry - as the field has come to be known - are just well-known theorems from other 
fields, but now in a different guise. This article just reinforces that opinion. 

2. Power diagrams and tropical hypersurfaces. 

In this section we show the equivalence of two definitions of the power diagram. First 
we discuss the original notion of Aurenhammer, then we show that these are equivalent 
to the definitions of Passare et. ai. 

2.1. The definition of Aurenhammer. The original definition of a power diagram 
is by means of distance function with weights. Take a point set {Pi, ■ ■ ■ , Pn} C. M". 
Throughout this article we assume that dim(CH({Pi, ■ ■ ■ , Pn})) = n. Assign a weight Wi 
to each point Pi. Then write down the functions: 

1 1 

(1) gi{x) = -\\x - Pi\\ - -Wi g (x) = min^ gi{x) 

Here we have used the following notation: 

n 
1=1 

Definition 1. For a subset a C {Pi, ■ ■ ■ , Pat} the set Pow(a) is the closure of 

{x G M" I gi{x) = g{x) Pi E a and gj{x) > g{x) Pj ^ a} 

The sets Pow({Pj}) are sometimes called cells. 

The separator Sep({Pj, Pj}) between Pi and Pj is the hyperplane. 

{xeW I giix) = gj{x)} 

It is no restriction to assume that Wi > 0. We may add some number to all of the 
functions gi'. the cells will not change. In fact we can add any function h: M" ^ M to the 
gi. We will get the same power diagram. Power diagrams are polyhedral subdivisions. 
We repeat definition 1 in |PR04j . 
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Definition 2. A polyhedral subdivision T of a polyhedron K C is a subdivision of 
K in polyhedra Ki, such that 

• The union of all sets Ki is K . 

• If Ki and Kj are both in T then so is their intersection. 

• Every compact subset L of K intersects only a finite number of the Ki. 

Theorem 1. The sets Pow(a) for a C {Pi, ■ ■ ■ , Pn} are a polyhedral subdivision o/R"". 
This polyhedral subdivision is called the power diagram. 

Proof. The cells cover M". They are the intersection of a finite number of halfspaces of 
the form: {x G \ gi{x) < gj{x)}. So they are polyhedra. □ 

Besides many similarities , there are a few crucial differences between power diagrams 
and Voronoi diagrams. In a power diagram the cell of a point may well be empty. It 
might happen that gi{x) > gj{x), for all i ^ j and all x G M". 

Also, whereas in Voronoi diagrams. Pi is always contained in its own cell, in a power 
diagram the cell of Pi might not be empty and still Pi does not lie in its own cell! 

The power diagram can be constructed much like the Voronoi diagram. Recall that 
the Voronoi diagram is the upper convex hull to the tangent planes to a parabola Xq = 
^iYll=i ^1)- -^^^ '^^^^ of power diagrams we do much the same thing. Look at figure 
HJ We take cylinders with radius around the lines x = Pi in M""*"^. The intersection 
of the cylinders with the parabola is a torsion zero space curve. We draw the planes in 
which they lie ( left figure ) and we consider the upper convex hull of these planes (right 
figure). The projection of its singular sides is the power diagram. 




Figure 1. Power lifting to the parabola 
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The above construction is also described in |Aur87j . section 4.1. When the = the 
construction reduces to the well-known construction of Voronoi diagrams. In that case, 
of course, one has to take the tangent planes to the paraboloid. 

2.2. Tropical hypersurface. The power diagram can be defined using only affine func- 
tions. With amoebas in mind, this was done in |PR04j . With tropical hypersurfaces in 
mind, this was done in |RGST03] . The separator between Pi and Pj is a hyperplane. It 
is given by 

x) = gj{x) ^ {x, Pi) = (x, Pj) 

We can define the power diagram solely using the affine functions fi{x) = {x,Pi) + Ci, 
where the coefficient Cj is 

(2) Q = 

Consider the function 

f{x) = max fi{x) 

i=l,... ,Af 

The following proposition shows that we can define the power diagram solely using affine 
functions. 

Proposition 1. For a subset a c {Pi, ■ ■ ■ , Pn} the set Pow(a) is the closure of 

{x I fi{x) = f{x) Pi e a and fj{x) < f{x) Pj ^ a} 

Remark 1. In |RGST03j . / is defined as the minimum of the functions fi. Taking the 
minimum or taking the maximum is just a convention though. 

The affine definition of power diagrams is more elegant. We lose the immediate appeal 
of the construction in figure ^ but two statements are now immediately clear. 

Proposition 2. The intersection of a hyperplane with a power diagram is again a power 
diagram. The image of a power diagram under an affine map M" is again a power 

diagram. 

So a power diagram can be sliced, after which we obtain a new power diagram. In fact, 
every power diagram in is a slice of a Voronoi diagram in M""*"^. 

Proposition 3. Let T he a power diagram in M". There is a Voronoi diagram T in 
M""^^ and a hyperplane H such that H (iT = T . 

Proof. We will explicitly construct such an T. Let {Pi, ■ ■ ■ , P/v} be the set of points in 
that determine T. Write 

Qi = {Pi,P^) x = {x,x') 

We have 

\\x-Qi\\ = Ik -Pill + {x' -Pi 
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So when x' = we get 

\\^-Q^\\ = \\X-P^\\+P^^ 

Now let the power diagram be given by functions Qi as in equation (^. It is no restriction 
to assume that Wi < 0, because only the differences Wi — Wj matter for the power diagram. 
Thus we can choose P/ = \/—Wi = ri. □ 

2.3. The Legendre transform. Dual to the Voronoi diagram is the Delaunay triangu- 
lation. In a similar way there exists a dual of the power diagram. This dual cell complex 
is not always a triangulation, similar to the case with the Delaunay triangulation. If we 
use the tropical point of view ( with affine functions ) we can best explain the dual object 
using the Legendre transform, see jHor94j . which we explain below. This dual object is 
also considered in j(TKZ94j . where it has the name coherent triangulation. 

We start with some definitions. Let P be a polyhedron in M""*"^. Let v be a vector in 
M". The lower faces of P with respect to v are those faces F of P such that 

Vx e P VA e M>o : X - At; ^ P 

A polyhedral subdivision of a polytope in M" is called regular if it is the projection of 
the lower faces of a polytope in M"^^. Not every polyhedral subdivision is regular, see 
chapter 5 in |Zie95j . or chapter 7 in |GKZ94j . 

Definition 3. The Legendre transform of a convex function f, with domain D cW^ is 

/(0 = sup((e,x)-/(x)) 

When the supremum does not exists, we put f{^) = oo. The domain Dom(/) of f are 
those ^ for which /(^) < oo. 

The Legendre transform of /(x) = |||x|| is the function itself. The Legendre transform 
of a linear function /j = (x, Pj) + Cj is < oo only when ^ = Pj. Theorem 2.2.7 of |Hor94j 
reads: 

Theorem 2. Let f = sup^g^ fa{x) be the maximum of a number of lower semi- continuous 
convex function. Then f is also a lower semi- continuous convex function. Furthermore 
f is the infimum over all finite sums: 



With theorem 121 we calculate the Legendre transform of the function that determines 
the power diagram. We have 

/ = max (x. Pi) + Ci 

l<i<N 

The domain where / < oo is the convex hull of the points {Pi, ■ ■ ■ , P/v}. We have 

/,(P,) = -Q^/(0=inf f-J^A, 



6 



D. SIERSMA AND M. VAN MANEN 



where the infimum is taken over all Aj such that 

N 

AjPj = ^ and Aj > and A,, = 1 

l<i<N 



It is no restriction to assume that only < n + 1 of the Aj are non-zero because any 
point in the convex hull of the Pi can be expressed as a sum of dimCH({Pi, ■ ■ ■ , P/v}) 
of the Pi. The infima are best thought of in a geometric way. Take the convex hull 
CH({(Pj, —Ci)}fLi). The lower faces form the graph of / over the convex hull CH({Pi, ■ ■ ■ , Pat}). 
We summarize our discussion in a theorem. 

Theorem 3. Let T be a power diagram. Let T be the lower convex hull of the lifted points 
[Pi, —Ci) wrt. to the vector (0, ■ ■ ■ , 0, 1). The polyhedral complex T is the graph of the Le- 
gendre transform f of f . The domain Dom(/) of f is the convex hull CH({Pi, ■ ■ ■ , P/v}). 

The geometrical construction is directly related to the construction of figure ^ If we 
put the upper convex hull of theorem |21 in figure [H we get figure |21 What exactly happens 




Figure 2. The power lift and the Legendre transform 

here can best understood by writing out the equations. In figure ^ the graph of xq = 

is drawn. Then the cylinders in the picture are ||x — Pj|| = rf. Or ||x — Pj|| = Wi. Hence 

the three planes of which we determine the upper hull are 

(3) xo = (x, Pj) + 

The triangle in figure |21 is the plane Xq = — X] -^jCj where the Aj are defined by x = ^ AjPj. 
In the point Pi we have that the graph of /(x) = maxj /(x) lies ^ above the paraboloid 
xo = whilst the graph of f\^) lies ^ below the paraboloid there. 
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Theorem |31 gives a dual triangulation of the power diagram. 

Definition 4. Given a power diagram T specified by N affine functions, the coherent 
triangulation a is the set of points ^ G CH({Pi, ■ ■ ■ ,-P/v}) where f is not differentiahle. 

Since the coherent triangulation is the projection of the lower convex hull, it has also 
the structure of a polyhedral complex. In the case when all the weights Wi are zero this 
coherent triangulation is of course the Delaunay triangulation. 

A number of points in {Pi, ■ ■ ■ , Pat} in is affinely in general position when the 
convex hull of each subset of length n + 1 of {Pi, ■ ■ ■ , P/v} is full dimensional. For Voronoi 
diagrams this definition is not sufficient. It might still be that n + 2 lie on a sphere. 
In that case the Delaunay triangulation is not a simplicial complex. For the Delaunay 
triangulation to really be a triangulation we need that the points {Qi, ■ ■ ■ , Qtv} are affinely 
general position in M""*"^ where 

Q, = {pJ-\\P,\\) 

Definition 5. A power diagram f = max/j is in general position when the points Qi 
are affinely in general position in M""*"'^. Here 

Qi = {Pi,-Ci), l<i<N 

The following theorem is left to the reader. 

Theorem 4. The coherent triangulation is a simplicial complex when the power diagram 
is in general position. 

Consider the function 

This function hides the Gateau differential. Namely, let /i : M" ^ R be a function. If the 
limit 

/z-(x;0 = lim^^^ + ^^^-^^^^ 

^ ^ m t 

exists, it is called the Gateau differential. See |Hor94j . theorem 2.1.22. For a convex set 
K C M" the supporting function of K is the function 

^ sup(x,0 

Theorem 2.2.11. in |Hor94| says that the Gateau differential ^ f'{x;C,) is the supporting 
function of 

{e I P(x,O = 0} 

Let us see an example of what that means. Take to be one of the points of {Pi, ■ ■ ■ , P/v}. 
Then /(Pj) = — Cj. So the set of a; where F{x,C,) = are those for which f{x) = {x, Pi)+Ci. 

The Gateau differential is called Clarke's generalized derivative in |APS97j . The set of 
which it is the support function is called df{x). It is stated in that article that df{x) is 
the convex hull of the gradients of the functions fi for which fi{x) = f{x). 
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That last statement and the theorem 2.1.22 are all equivalent to what is neatly formu- 
lated in proposition 1 of |PR04j : 

Theorem 5. There is a subdivision of the convex hull CH({Pi, ■ ■ ■ ,P/v}); dual to the 
power diagram T . The cell Del(a) dual to Pow(a) & T is 

Del(a) = I F{x, = Vx G Pow(«)} 

Reversely 

Pow(a) = {x I = e Del(a)} 

Remark 2. Similar statements are in the review by Pennaneac'h, see |Pen03j . However 
Pennaneac'h does not discuss the relation between the figures U\ and{^ Also she does 
not make clear what the relationship between this theorem and the well-known theory of 
Voronoi diagrams/ Delaunay tesselations is. Her review however is heartily recommended. 



3. Disappearing faces 

We have noted before that faces in the power diagram might disappear. This will 
happen for some Pj G {Pi, ■ ■ ■ , P/v} iff. 

fj{x) < f{x) Vx G W 

Cells in the power diagram correspond to vertices in the coherent triangulation. In that 
case it is very easy to see which vertices disappear. They are the ones lying in the epigraph 
of the Legendre transform. Remember that the epigraph of a function h: D ^ M is defined 
as: 

I h{x) <t X e D} 

Let us sum up some consequences: 

If Pj lies outside the convex hull of CH({Pi, ■ ■ ■ , P^} \ {P,}) then the cell Pow({P,}) is 
not empty. Proof: the lifted point [Pj, —Cj) forcibly is part of the lower hull of the points 

The transition on the power diagrams where a cell disappears is exemplified in figure |21 
In the pictures the convex hull of the points (Pj, — Cj) is drawn. The graph of the Legendre 
transform / is formed by the lower faces of that polytope. In the ground plane we see the 
projection of those lower faces: the coherent triangulation dual to the power diagram. In 
the ground plane also the circles with radius ^ywi are drawn. 

On the left hand side of figure El we see a point with a small circle around it. This point 
disappears when we make one of the circles in the pictures ( the one on the upper right 
) bigger. Such has been done on the left hand side of figure ^2 The upper right circle 
is much bigger there and the lifted point corresponding to the small circle lies inside the 
epigraph of the Legendre transform /. We thus see that there are three sorts of events on 
planar coherent triangulations: 

• Convex hull events and edge flips 

• Disappearing vertices. 
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Figure 3. Transition on power diagrams 



Convex hull events and edge flips already occur in families of Voronoi diagrams. A detailed 
study of these phenomena can be found in |Lin02j . However in a dual triangulation of a 
Voronoi diagram, none of the cells Pow({Pj}) can be empty. 



4. The k-TU order power diagram and hyperplane arrangements 

Given a power diagram we define a kth order power diagram. Consider all subsets I of 
length k in {!,■■■ , N}. Denote 



i&I 



and 



Denote also 



So /(I) = /. 



/(''^ = max // 

/C{l,--,Af}|/|=fc 



Definition 6. The k-th order power diagram Pow*-^-* is the power diagram of the fj. 
Denote Del'-'^-' its dual. 

Our definition is admittedly rather arbitrary. One can take any weighted sum of the 
points. 

Lemma 1. The k-th order power diagram is a polyhedral complex. The point x lies in the 
cell I of the k-th power diagram iff. fi{x) > fj{x) \li E I and Vj G {1, ■ ■ ■ , A^} \ /. 
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Proof. The first part of the lemma is obvious. The second part of the lemma boils down 
to the following obvious statement. Let {ai}^i be a set of real numbers. If / is the index 
set of the k biggest of the {ai}^i then aj is the biggest of the numbers 

{ij^a, I \J\ = k} 

□ 

There are two sources of cells in Pow'-'^'^^^. 

• Let Pow('=)(/i) and Pow(^')(/2) be two power cells in the A;-th order power diagram. 
Assume that the two cells are adjacent, that is there is an edge between the two 
vertices in the dual coherent triangulation Del*^^^. This means that /i n/2 has k + 1 
distinct elements. Thus the cell of /i fl I2 in Pow*^*^"^^^ is not empty. 

• Let the cell of / C {Pi, ■ ■ ■ , P/v} be an empty cell. That is: 

fi{x)<fj{x) ^xeW WJc{Pu--- ,Pn} \J\ = k 

but that at some point Xq G M" there is only one j G {Pi, ■ ■ ■ , P/v} \ / for which 
fj{xo) is bigger than all /i(xo) for i E L 




Figure 4. A sample of three different kinds of Del Del*-^-*. If there is 
a disappeared vertex in Del that vertex may show up in Del^^). The dots 
indicate such vertices. 



Theorem 6. The union of the power diagrams of order k is a hyperplane arrangement. 
The hyperplanes are the separators. Each such separator is equipped with a polyhedral 
subdivision, that is: the separator between Pj and Pj is a hyperplane V divided in at most 
N closed convex subsets Ki. If x lies in the interior of Ki then x lies on the l-th order 
power diagram. 

Proof. That there is a polyhedral subdivision on the separator is rather obvious. The 
points of the /c-th diagram is the intersection of two closed convex subsets, so its intersec- 
tion is also convex. □ 

The aforementioned hyperplane arrangement is not a generic arrangement. Generically, 
at most n hyperplanes intersect. In this case, n + 1 hyperplanes intersect when they 
correspond to some cell in a higher order diagram. A hyperplane arrangement does not 
determine in general a power diagram. 

Drawings that look like a power diagram need not be one, just as not all triangulations 
are coherent. For examples of the latter phenomenon see for instance section 5.1 in |Zie95j . 



where coherent triangulations are called regular. As an example of the former, see figure 
m Not all directions will lead to an intersection of the three arrows inside the cell. 
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Figure IHFour regions, marked by Pi, P2 , 
P3 and P4. Qi = (1,-2), Q2 = (2,2) and 
Q3 = (—2, 0). If the halflines from the Qi are 
allowed to move in the specified directions 
we always have a polyhedral subdivision, but 
only a power diagram when the back ends of 
the arrows intersect inside the triangle. 



5. The Morse poset 

In jSvM04j we introduced the Morse poset for Voronoi diagrams. Here we generalize 
the Morse poset, so that it is defined for power diagrams. 
Consider the function 

q = min Qiix) 

This function is a topological Morse function. When does it have a critical point? 

Theorem 7. The function g does not have critical points outside the convex hull of 
{Pi,--- ,Pn}- 

Proof. Let x be a point outside the convex hull of {Pi, ■ ■ ■ ,P/v}. Hyperplanes in T^M"' 
can be identified with affine hyperplanes in M" that pass through x. Choose a hyperplane 
H in X that is parallel to some hyperplane that separates x from CH({Pi, ■ ■ ■ , P/v}). It 
follows that all gradients ^ = x — Pi point outward from x at one side of H. Thus the 
function miuj Qi does not have a critical point at x. □ 

Consider the coherent triangulation dual to a power diagram. In the special case of a 
Delaunay triangulation we have defined a poset of active simplices in the triangulation. 
They were defined using the critical points of mingfj. We repeat the procedure for power 
diagrams. 

First define the separator of a cell a G Del({Pi, ■ ■ ■ ,P/v}). The separator Sep(a) of 
a is the intersection of all the separators of Pj and Pj such that {Pj,Pj} C a. If a is a 
vertex then the separator is defined to be just the whole of M". Also, for a polyhedron 
P in W^'^^ the relative interior is the interior of the polyhedron as a subset of the affine 
span Aff(P) of P. If dim a > then put 

(4) c{a) = Sep(a) n Afr(a) 

For a vertex put c{{i}) = Pi. We call the simplex in the coherent triangulation active if 

• c{a) is contained in the relative interior Interior(CII(a)) of Cll(a), and 

• g{c{a)) = gi{c{a)) for some Pi E a 

With this definition it might well happen that a vertex is not active. Here is a typical 
counter intuitive example. Take three points, say Pi = (—1,0), P2 = (1,1) and P3 = 




12 



D. SIERSMA AND M. VAN MANEN 



(1,-1). Then place a large circle around Pi: ri = 5. Place two smaller circles around 
P2 and P3: r2 = and = y/|. If we look from below at the graph of mingfj we see 
a picture alike the one in figure El The Morse poset of {(Pj, u'j)}j=i,... ,7v consists of the 



01 
-50l 
-lOOl 



Figure 5. Typical counterintuitive example 

active subsets of the coherent triangulation. We write 

ga{x) = mm gi{x) and = max fi{x) 



Theorem 8. There is a one to one correspondence between critical points o/miuj^fj and 
active cells in the coherent triangulation. An active cell of dimension d corresponds to a 
critical point of index d. 

Proof. The function g{x) = mmigi{x) can only have a critical point at xq when at least 
two of the (7i(xo) are equal to g{x). The index set for which gi{xo) = g{xo) is the face a. 
The function can also only have a critical point when xq G CH(a). We can thus assume 
Xq = c{a). 

The index of the critical point can be calculated with theorem 2.3 in |APS97j . First 
note that for all x G M" the Hessian of gi{x) is always the unit matrix I. 
Consider the subspace 

T{x) = Hig^ker "^(^^o) = njg„ker(a; - Pi) 

Clearly T(c(a)) has the same dimension as Sep(a). The dimension of T is n — dim a. 

The index of the critical point is n minus the dimension of T(x), it is thus the dimension 
of a. 

Conversely if a is active then ga restricted to CH(a) has a supremum at c(a). 

g{x) < g^x) < g4c{a)) Wx G CB.{a) 
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Also there is a neighborhood U of c{a) where g{x) = 5'a(c(a)) and so g has a critical point 
of index dim(Q!). 

Note that the reasoning is independent of genericity conditions. □ 

Finally we remark that though power diagrams are affinely defined, activity is not 
retained under affine volume preserving linear transformations. Through such transfor- 
mations the power diagram and its dual coherent triangulation do not change. The Morse 
poset though, can change drastically after an affine transformation. 



6. Discrete Morse theory 

We recall the main definitions of the discrete Morse theory developed by Forman. A 
very good introduction to the subject is |For02j . We show that, though the coherent 
triangulation is not in general a simplicial complex, it has another nice property: for 
every two faces a C 7 with 2 + dim(a) = dim(7) all there are two faces f3 of the coher- 
ent triangulation in between a and 7. Hence, discrete Morse functions on the coherent 
triangulation have the same nice properties as discrete Morse functions on a simplicial 
complex. 

Let T be a polyhedral subdivision. By orientating all the faces T becomes a complex. 
Definition 7. A function h: T —>■ M. is called a discrete Morse function if for all P E T 



(5) #{a G r I 1 + dim(a) = dim(/?) aGp h{a) > h{(3)} < 1 
and 

(6) #{« e T I dim(a) = 1 + dim(/5) p Ca h{a) < h{p)} < 1 
In case both numbers are zero for some (3 & T , (3 is called critical. 

Analogous to the case of a power diagram and its dual coherent triangulation, we can 
define a Morse poset of a discrete Morse function h as the set of critical faces of h. 

If the power diagram is in general position T is a simplicial complex. Discrete Morse 
function on a simplicial complex have a very nice property, see lemma 2.6 in |For98j . If 
a is a face of (3 and /5 is a face of 7. Label the vertices: 



a 



{Pi,---,P,.} /3 = {Pi,--- ,P,+i} 7 = {A,---,Pfe+2} 



Because we are dealing with a simplicial complex there is another face between a and 7: 
^ = {Pi, ■ ■ ■ , Pfc, Pfc+2}: 



a 



I5 





{Pi^ ■ ■ ■ 5 Pk+2] 



Pk. Pi 
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Let / be a discrete Morse function on the simplicial complex. Suppose that /(«) > 
fiP) > fil), that is both ® and © hold. We also have f{a) < f{5) < f{-f), or 

/(7) > f{S) > f{a) > f{(3) > /(7) 

This is a contradiction, so we see that on simplicial complexes at most one (0) and © 
can hold. 

However this argument is not restricted to simplicial complexes. All that is needed is 
that if a and 7 are two cells with a C 7 and 2 + dim a = dim 7 then there are at least 
two cells in between a and 7. 

Theorem 9. For a coherent triangulation a as in the above between any two cells a and 
7 with 2 + dim a = dim 7 and a C 7 there are j3i and (32, both cells in the coherent 
triangulation, both having dimension 1 + dim a such that a (Z jSt (Z j. 

Proof. Consider the cone Cone(7, a) of 7 over a: 

Cone(7, a) = {(3 \ a Z f3 C 7} 

Intersect that cone with the n — k + 1 dimensional plane through zero and orthogonal to 
the k — 1-dimensional plane affinely spanned by a. The intersection is a cone in the two 
dimensional plane. It has two extremal vectors .^1 and .^2- These correspond to the faces 
f3i and we were looking for. □ 

Hence even if the coherent triangulation we have is not a triangulation then we still 
have that not both (0) and (jH)) can hold. 

If we mod out Cone(7, a) by the directions Aff (a) along a we get a vertex of a polyhe- 
dron and the faces between a and 7 are the edges of a connected fan. So if we have that 
there can be any number of and /5j2 between a and 7. 

r/?ii c An 

« ^ < : : : > 7 

[Pii c AiJ 

On any coherent triangulation there exists a discrete Morse function: h: a ^ dim a. 
So let us have a discrete Morse function h on the coherent triangulation. The gradient 
of h consists of arrows « — > /3, drawn when a G f3 and h{a) > h{(3). In general such 
a collection of arrows is called a discrete vectorfield. A V^-path is a sequence of arrows 
such that the simplex pointed to contains a second simplex that points to another one, of 
higher dimension. A closed V-path is a circular \/-path as in equation (|7j). Theorem 3.5 
in |For02j states that a discrete vectorfield is the gradient of a discrete Morse function if 
and only if it contains no closed ^-paths. 

(7) ai ^ /?i D 02 ^ /52 D ■ ■ ■ D Ofc ^ D ai 

The dual to the coherent triangulation, the power diagram, is also a polyhedral complex. 
Given a discrete Morse function g on the coherent triangulation, clearly — g is a discrete 
Morse function on the power diagram. 
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7 7 




Figure 6. A cone of simplices filled up with a discrete vectorfield 

7. From Morse poset to a discrete vectorfield 

In this section we will prove the main new theorem of this article. 

Theorem 10. On Del({Pi, ■ ■ ■ ,-P/v}) there exists a discrete Morse function h such that 
the Morse poset of h equals the Morse poset of g{x) = mini<j<7v 5fj(x). 

In fact we will not construct h directly. Instead we will construct a discrete vectorfield 
without closed l^-paths. In the reasoning below we disregard the disappearing faces. 
A face a can be non-active for one of two reasons 

Down The separator Sep(a) can lie outside Interior(CH(a)), or if that is not the case: 
Up Some point P.^ G {Pi, ■ ■ ■ , Pn} \ ol can lie closer to c(a) then the points in a, i.e. 

(8) gi{c{a)) < gj{c{(x)) Vj G a 

Suppose that a is not active for the "Down" reason. Then in the case of Voronoi diagrams 
a cannot be a vertex or an edge. In general with power diagrams a cannot be a vertex. 

Suppose that a is not active for the "Up" reason and that dimDel(a) = n. Then we 
have 

{c(a)} = n,e,Pow({P,}): 
we cannot have (jH}. We conclude that no a, not active for the "Up"-reason has dim Del(a) = 
n. 

We start by showing that for every face, not active for the "Down" reason, there is a 
unique lower dimensional simplex not active for the "Up" reason. Then we construct for 
each face not active for the "Up" reason a part of the discrete vectorfield. 

Lemma 2. Let f3 G Del({Pi, ■ ■ ■ , Pat}) be a face not active for the "Down" -reason. Then 
the closest point to CH(/3) from c(/3) is c{a) for some proper face a of p. At c{a) we have 

(9) g,{c{a)) < gMc^)) Vj G /3 \ « 

Hence a is not active for the "Up" -reason. Conversely if holds for j3 D a then (3 is 
not active for the "Down" -reason. 

Proof. We start with "^" . 

A special case we first have to handle is when c{(3) lies on the relative boundary of 
CH(/5). This means that c(/3) G CH(a) for some proper face a of p. Because Sep(/3) C 
Sep(Q;) it follows that c(/3) = c{a). Equation (jH)) follows automatically. 



16 



D. SIERSMA AND M. VAN MANEN 



Now we can safely assume that the distance from c(/3) to CH(/3) is > 0. Denote y the 
closest point on CH(/5) from c{(3). If a is a vertex we are done, so suppose that a is not 
a vertex. The line from y to c{j3) is orthogonal to Aff(a) and hence parallel to Sep(a). 
Clearly c(/9) lies in Sep(a). So y also lies in Sep(a). But y also lies in Aff(a). So, by the 
definition of c{a) in equation (jH) we have that y = c(a). 

Because CH(/3) is a convex set the hyperplane with normal c(/3) —c{a) passing through 
c(a) separates Interior(Cif(/5)) from c(/3). Thus the angle to c{a) to c(/3) must be 
obtuse. Consequently: 

(10) \\c{/3) - PkW > ||c(/3) - c{a)W + ||c(a) - P^W 
if Pk is a vertex in (3 D a. By Pythagoras, equality holds if G a: 

(11) l|c(/5) - PkW = l|c(/3) - c(«)|| + ||c(a) - PkW 

Then equation (jiup becomes, taking the factor ^ we used in the definition of g in equation 
(HI) into account: 

(12) ^,(c(/5))>i||c(/5)-c(a)|| + (7fe(c(«)) ke(3\a 

Again, equality holds when k E a. Putting this together we get Q, and we see that a is 
not active for the "Up" -reason. 

Next do the part. We assume Q and a G (3. From we get for j G /3 \ a: 

(13) gM»)) + ^||c(/3) - c{a)W < gM»)) + ^Wm - c(a)|| = gMP)) = 9MP)) 

The special case to handle first is where equality holds in (jl3|) . Then c{/3) = c{a) and so 
c{f3) is not an element of the relative interior Interior(CH(/5)). So f3 is not active for the 
"Down" -reason. 

Now we can safely assume strict inequality in (fT!?|l . The simplest case is a /5 for which 
f3 = aU{j}. The three points Pj, c{f3) and c{a) all lie in AS{(3). From (O we get (HU)) 
with strict inequality. The line segment from c(a) to c{(3) is thus an outward normal to 
CH(/5) in Aff(/5). Consequently is not active for the "Down" -reason. The general case 
is similar. □ 

Let again a be not active for the "Up" -reason. Denote K the set of indices for which 
(jHI) holds. Denote Up(a) the set of simplices: 

(14) Up(a) = {/3 G Del({Pi, ■ ■ ■ , Pn}) \ aC(3caUK} 

Lemma 121 characterizes the elements of Up(a) as those /3 D a for which c{a) is the closest 
point on CH(/?) from c(/3). The complement of the Morse poset is divided into different 
subsets Up(a), one for each face not active for the "Up"-reason. 

In Up(a) identify those faces 6 such that there are no 7 in Up(a) with 6 a proper face of 
7. Then we see that Up(a) contains different cones Cone{6,a). These cones can be filled 
up by a discrete vectorfields, as is shown in figure El It is possible to take a number of 
simplicial cones such that no discrete vectorfield exists that fills up these cones completely. 
For instance, take two triangles: {Pi,P2}, {Pi, P2, P3} and {Pi, P2, P4}. Obviously there 
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is no way to fill these cones up with a discrete vectorfield. Therefore, we need to prove 
that Up(a) can be filled up with a discrete vectorfield. 

We will establish another criterion for which a /5 G Del({Pi,- ■■ ,Pn}) with a C /5 
is an element of Up(a). Take the point c(a) outside the polytope Pow(a) in Sep(a) 
and consider the faces of Pow(a) for which the outward pointing normal makes an angle 
smaller than 90 degrees with x — c{a). Then we get a number of faces ■ ■ ■ , Pi in 
the Delaunay triangulation that we see as faces Pow{j3i), ■ ■ ■ ,Pow(/3;) on the relative 
boundary 9(Pow(a)) of Pow(ci;). This is illustrated in figure [7| A precise formulation is 




Figure 7. View from c(a) to Pow(a) 
contained in the following lemma. 

Lemma 3. Let (3 D a be an element o/Del({Pi, ■ ■ ■ , Pn})- Let x be a point in the interior 
o/Pow(/5). //, either X = c{a) , or 

(1) the line through x and c{a) intersects Pow(a) in a segment, and 

(2) the segment [x, c(a)] contains no point of the interior o/Pow(a) 

then (3 G Up(a). Conversely, if 1 and 2 hold for some P D a, with a not active for the 
"Up" -reason, then /3 G Up(a). 

Proof. We first need to treat the case x = c{a). In that case c(a) lies on the relative 
boundary of Pow(a). So x G Pow(/3), P = aU L And fj{c{a)) = fa{c{a)), for all j G /. 
Obviously Up(a;) consists exactly of all those 7 G Del({Pi, ■ ■ ■ , Pn}) with 

a C'j C p = aU I 

Now we can assume x 7^ c{a). Again x lies on the relative boundary of Pow(/?) and 
/3 = a U / for some maximal set /. We put xt = tx + {1 — t)c{a). Hence xi = x and 
xo = c(a). It also follows for j G / 

f,{xt) = tf,{x) + (1 - t)/,(c(a)) = tfM + (1 - t)fM»)) 

fA^t) - faiXt) = (1 - t)(/,(c(«)) - /«(C(«))) 

Now we have 

fj{c{a)) > faic{a)) ^ 3t > 1 : fj{xt) < faixt) 

So that 

3t> 1 XtE Pow(a) fjic{a)) > fa{c{a)) 
And this is exactly what we needed to prove. 

Now suppose that P D a, and a not active for the "Up" -reason, but that 1 and 2 do 
not hold. 
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The only point of Pow(q;) on the hne through x and c{a) is x. It follows that (3 \a 
contains more than one index, i.e. /3 = a U {ji, ■ ■ ■ , >} with r > 2. In a sufficiently 
small neighborhood of x the halfspaces Hi = {fj. < fa} define Pow(a) as a polyhedron 
in Sep(Q;). Because {xt | t G M} touches Pow(a) there exist an index k, say k = ji, and 
a t < 1 such that fj{xt) — fa{xt) and hence fk{c{a)) < fa{c{a)). And so /5 D a is not an 
element of Up (a). □ 

From lemma 121 we see that we can fill up Up(a) with arrows of a discrete vectorfield. 
Indeed, with its dual on Pow(a) we can obviously do this because it is contractible, so we 
can also do it with Up (a). 

We conclude that the complement of the Morse poset can be divided into separate parts 
Up(a), for all a not active for the "Up" -reason. Each of these parts cones represents what 
Matveev calls a polyhedral collapse in |Mat03j . To show that we indeed get a discrete 
vectorfield we need to prove one more thing: there should be no closed ^-path. 

Suppose that we have a closed ^-path. Inside one Up(Q!) there is no such closed V^-path, 
so at some simplex the V^-path should jump from one Up(7i) to another Up(72), as in 
figure IHl 




7i = "1 



Figure 8. A V-path, with a jump from Up(7i) to Up(72). 



Lemma 4. Suppose that we have a V-path that jumps from Up(7i) to Up(72) as in equation 

m- 

ai ^ Pi C Cone(7i U /i,7i) 

(15) U 

ai+i (3i+i C Cone(72 U 12,72) 

Then ^^,(0(71)) > ^^3(0(72)). 

Proof. Because 71 is a proper subset of 71 U 72 we see that: 0(71 U 72) ^ Interior (CH (71 U 
72)). The closest point to CH(7i U72) from 0(71 U72) is 0(71), as we have shown in lemma 
121 In particular 

(16) ||c(7i U 72) - c(7i)|| < ||c(7i U72) - 0(72)11 

Here, the inequality is strict. The closest point is unique, because CH(7i U 72) is convex. 
For the reasons explained in the proof of lemma |21 we have 

(17) ^^,u72(c(7i U72)) = ^7i(c(7i U72)) = ^l|c(7iU72) -c(7i)|| +^^,(0(71)) 
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The same identity holds with 71 and 72 exchanged, thus: 

(18) ^l|c(7iU72) -c(7i)|| +^7^,(c(7i)) = ^||c(7i U 72) - c(72)l| + ^772^72)) 
Putting (fTHIl and |TH|) together we get the desired resuh. □ 

Lemma 0] says that to each Up(a) appearing in the \^-path we can associate a number, 
and passing from one Up(a) to another that number strictly decreases. Hence we can not 
have a closed l^-path. The proof of theorem El is complete. 

Remark 3. We are not entirely satisfied with the proof of our main theorem. What is 
lost in the proof is a direct relation between g and the discrete Morse function that gives 
the discrete vectorfield we construct. We tried using a ^'^Pxecu{a) 9{^) ^ function 
on the coherent triangulation hut this approach did not give the desired results. 

Remark 4. We were neither entirely satisfied with discrete Morse theory, because it seems 
to he ill equipped for dealing with polyhedral collapses: our vectorfield is not unique. The 
non-uniqueness only lies in the choice of the arrows possible in Up(a). One could also 
try to develop a generalization of discrete Morse theory in order to include polyhedral 
collapses. 

8. Examples 

Let us return to the counterintuitive example of figure El There the Morse poset is as 
in figure IHl In the picture we see three cones, consisting each of two faces in the simplest 
triangulation one can think of. 




Figure 9. Morse poset and discrete vectorfield of figure El 

Another example is the one called 4300/ in our previous article |SvM04j . Let the active 
subsets be {{Pi, P2}, {A, ^^3}, {A, A}}- Then if a = {Pi, P2, ^"3, ^"4} the face {Pi,P4} 
is not active for the "Up"-reason. In that case both {Pi,P3,P4} and {Pi,P2,P4} and 
{Pi, P2, P3, P4} are not active for the "Down" reason. There is one cone consisting of the 
polyhedral collapse: 

{Pi, P4} ^ {Pi, P2, P4} {Pi, P3, P4} ^ {Pi, P2, P3, Pa} 

The decomposition of that cone into arrows is not unique. We might just as well write 

{Pi, P4} ^ {Pi, P3, P4} {Pi, P2, P4} - {Pi, P2, P3, Pa} 
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Another application is the classification of all Morse posets of Voronoi diagrams for 4 
points in the plane in |Sie99j . Unfortunately there were two cases missing, which are now 
included in figure ITUl 












Figure 10. 3 and 4 points in M^: all Morse posets from Voronoi diagrams 
with the accompanying discrete vectorfields 



9. The medial axis and Maslov dequantization 

It is well known that the medial axis of a compact embedded submanifold M in can 
be calculated as an approximation of a Voronoi diagram, see for instance |ABK98j . We 
use tropical geometry to give a new interpretation of that fact. Recall that the corner 
locus of a convex function / is the set of points where / is not different iable. 

We have seen in the above that the Voronoi diagram can also be calculated using the 
affine functions /j, instead of the distance function. The corner locus of the function 



f{x) = max(a;. Pi) 



is the Voronoi diagram. It is natural to try the same procedure with the embedding 
7: M — s> M". We will show next that the medial axis is a corner locus as well. 



Proposition 4. Let M be a compact smooth manifold, embedded by M 
Write for s e M: 



fsix) = (x,7(s)) - -his) 



The function f defined by 



(19) 



fix] 



def 



sup/s(x) 



lim 

ft— >oo 



log U,,hf^(^)ds) 
log{h) 



is a convex function, and satisfies the stated equality. Its corner locus is the medial axis 
ofM. 
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Proof. The function /(x) is convex, because it is the supremum of a number of convex 
functions. Because M is compact the supremum is attained for one or more points s = s{x) 
in M. We see that x G M"' and s e M , for which f{x) = fs{x), are related by 

(20) (x-7(s),7(s)) = 

I.e. X lies on the normal of one or more s G M. Of those that satisfy (PUI) we choose one 
that has the highest value of fs^ [x) , and thus the lowest value of 

1 1 

The s we get from sup^gj^,^ fg corresponds to the point closest to x on the submanifold 
M. And the points where /(x) is not differentiable form exactly the medial axis of M. It 
remains to prove the formula for sup^gj^^ fg in (fT^. 

For this we use the so-called Maslov dequantization. Maslov dequantization is a fancy 
name for the identity: 

(21) lim log;,(/i'' + /i*) = max(a, b) 

h—>oo 

that holds for any two real numbers a and b. □ 

For practical purposes this procedure works well. The integrals can in most cases not be 
calculated explicitly, but if one uses numerical integration one can plot /, as is illustrated 
in figure HH 




Figure 11. The graph of sup^gjj /s(x) when 7 is the ellipse. 



The Maslov dequantization as in equation (|2ip is strongly related to a construction in 
toric geometry, described in |Ful93j . section 4.2. Namely, let again 



f,{x) = {x,Pi) + Ci z = l,---,N 
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be affine functions defining a power diagram. Define the following function 

H:W-^W H(x) = ^^=J ' 

It is proved in |Ful93j that is a real analytic diffeomorphism from to the interior of 
CH({Pi, ■ ■ ■ , -Pat}), when that interior is not empty. One recognizes if as a gradient: 

N 



Now consider the gradient: 




i=l 



i=l 



Proposition 5. Suppose that the interior o/CH({Pi, ■ ■ ■ ,-P/v}) is not empty. The map 
Hh is a real analytic diffeomorphism from to the interior of CH({Pi, ■ ■ ■ , Pat})- 

Proof. The fi are of the form {x, Pi) +Cj. Put fi{u) = {u, Pi) + log{h)ci. We already know 
that 

u ^ H{u) = 

is a real analytic diffeomorphism from to the interior of CH({Pi, ■ ■ ■ , Pat}). Because 
we have Hh{x) = H{x\og{h)) the proof is complete. □ 

Taking infinitely many points and a Riemann sum we see that the gradient of the 
function log^ (^Jj^ h^'^^^ds) is a diffeomorphism from to the region in M" bounded by 
the compact embedded manifold M C M". Now in the limit h —* oo, the image of M 
is the medial axis. This is another way of saying that the dual of the medial axis in the 
sense of the theorem of Passare and Rullgard is the manifold M itself. 
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